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Abstract 

We conduct depth-resolved three-dimensional Direct Numerical Sim- 
ulations (DNS) of bi-disperse turbidity currents interacting with com- 
plex bottom topography in the form of a Gaussian bump. Several 
flow characteristics such as suspended particle mass, instantaneous wall 
shear stress, transient deposit height are shown via videos. Further- 
more, we investigate the influence of the obstacle on the vortical struc- 
ture and sedimentation of particles by comparing the results against 
the same setup and but with a flat bottom surface. We observe that 
the obstacle influences the deposition of coarse particles mainly in the 
vicinity of the obstacle due to lateral deflection, whereas for the sedi- 
mentation of fine particles the effects of topographical features are felt 
further downstream. The results shown in this fluid dynamics video 
help us develop a fundamental understanding of the dynamics of tur- 
bidity currents interacting with complex seafloor topography. 

1 Introduction 

Turbidity currents represent a large-scale geophysical flow phenomenon in 
which the suspension travels along the bottom surface due to a horizontal 
hydrostatic pressure gradient ([I]). These currents play an important role 
within the global sediment cycle, and in the formation of deep-sea hydro- 
carbon reservoirs. They can travel up to O (1, 000 km) along the bottom of 
the ocean and transport 0(km 3 ) of sediment into deep water regions. 
Due to the depositional/erosional nature of turbidity currents, the inter- 
action of such currents with seafloor topography can produce pronounced 



Figure 1: Schematic of the initial set-up of the lock-exchange configuration 
interacting with bottom topography in the form of a Gaussian bump. The 
full lock has dimensions L s xHxW = 1x2x3 and contains a bi-disperse 
suspension of density of p\. The obstacle has a height of h = H/A and is 
centered at (x^, Zb) = (5.5L S , W/2). At time zero, the membrane is removed 
and a turbidity current starts propagating along the bottom surface. 



features such as channels, levees, sediment waves and gullies. 
In this study, we present results from DNS calculations of a bi-disperse 
turbidity current interacting with a complex topography in the form of a 
Gaussian bump. The results are obtained with our in-house code TURBINS 
([2]), which allows us to perform large-scale simulations of these currents. 
We solve the Navier-Stokes and transport equations to describe the fluid 
motion and evolution of particle concentration, respectively. 



2 Results 

The suspension includes two particle sizes (see figure [I]), i.e. coarse and fine, 
with dimensionless settling velocities 

u c s = 0.03, (1) 
u f = 0.006, (2) 

and initial mass loading of 50% and 50%, respectively. The Reynolds num- 
ber, defined by 

Re=^ (3) 

is set to Re = 2,000. In equation O), the buoyancy velocity u?, serves as 
the characteristic velocity (note: "refers to any dimensional quantity). We 
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remark that in all results, the particle concentration c refers to the volumet- 
ric particle concentration divided by the initial particle concentration in the 
suspension in the lock region (see figure [TJ . 

The video presented [link here] consists of mainly two parts. The first part 
presents various features of the current produced by the lock-exchange con- 
figuration interacting with the complex topography. Visualization of various 
flow features such as particle concentration, particle suspended mass, deposit 
height, and wall shear stress help us gain better insight into the flow dy- 
namics and turbidity current characteristics. 

We further investigate the fate of both particle sizes by tracking 48, 000 La- 
grangian markers for each particle size colored by their initial location in 
the suspension in the lock region (see figure [TJ . From these results we ob- 
serve uniform mixing and deposition of fine particles as a result of current's 
turbulent motion. 

The second section focuses on understanding the influence of the complex 
topography on the current. Toward this goal, we compare our results ob- 
tained for the case with mounted obstacle against the case without any 
topographical feature, i.e., a flat bottom surface. 

As a result of the interaction of the current with the bump, we observe 
non- uniform vortical structures including separation regions, bi-section and 
re- attachment, and turbulent mixing as opposed to uniformly spaced longi- 
tudinal streaks (lobe-and-cleft instabilities) present in the flat surface simu- 
lation. 

Moreover, comparison of the final deposit profiles for both particle sizes sug- 
gests that coarse particles are increasingly deposited in the vicinity of the 
bump due to lateral deflection. For the fine particle sizes, however, we notice 
more pronounced differences further downstream of the Gaussian bump. 
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